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Abstract 

\ We show how one can implement any local quantum gate on specific qubits in an array of qubits 

o ■ 

^Sj ' by carrying adiabatically a Hamiltonian around a closed loop. We find the exact form of the loop 
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> 



o 
o 

:^ 



and the Hamiltonian for implementing general one and two qubits gates. Our method is analytical 
and is not based on numerical search in the space of all loops. 
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I. INTRODUCTION 



and its further development in 

mm 

concrete realization of holonomic quantum gates 
has attracted a lot of attention. The attraction comes mainly from the fact that this type of 

S.antum computation may be inherently stable against some local errors. The basic idea of 
is as follows. One takes an n dimensional Hamiltonian with a degenerate k dimensional 
eigenspace 

H{t)m) = e{t)m), (1) 

where {\i{t)), i = 1 k} are the instantenous degenerate eigenstates of the Hamiltonian. 
Carrying this hamiltonian adiabatically around a loop in the parameter space, induces a 
generalization of the Berry phase in the form of a unitary operator in this subspace. This 
unitary operator is given by 

U = Pe-^c^, (2) 
where P means path-ordered exponential, A is the connection one form 

{i{t)\A\3{t)):={i{t)\d\j{t)), t,j = l to k, (3) 

and C is the loop around which the Hamiltonian is moved in the parameter space. Here we 
are assuming that the basis vectors of the Hilbert space have been arranged so that the first 
k states span the degenerate eigenspace of the Hamiltonian. 

Of particular interest to us is a scenario first proposed from a geometrical point of view in 
The physical description of the idea of p] is as follows: To make the connection one 
form time independent, one takes a family of iso-spectral Hamiltonians on an n dimensional 
Hilbert space defined as follows: 

H{t) = e*^i/oe-*^, (4) 

where Hq = €j2i=i = ^Pqi and X is an anti-hermitian operator acting on the Hilbert 
space. Thus Pq is the projection operator for the k— dimensional eigenspace, and the 
operator X determines the adjoint action and hence the curve in the parameter space. Note 
that the curve has been parameterized by the real parameter f G [0,1] and if one requires 
that the curve be closed in the parameter space, one should have 

Po = e^Poe-^. (5) 



i in fl( 



Since the suggestion of Zanardi and Rasetti in [1'] on holonomic quantum computation 



It is always possible to choose the basis of the Hilbert space so that 

(6) 



V 0/ 



For the Hamiltonian (HI) we have: 



|z(t)) = e*^|z), (7) 
from which we obtain the time independent connection 

m\A\m) = {i\x\j), i,j = i ■■■ k. (8) 

Note that this equation by no means implies the A = X. It only means that the restriction 

of the operator X to the degenerate subspace is equal to A. 

Since X is constant, calculation of the path integral is quite simple and gives 

U = Pe~fc^ = e-^. (9) 

Note that the operator X determines the curve in the parameter space completely. However 
only its projection to the degenerate subspace is determined by the operator U = e^^ which 
we want to implement w holonomically through equation and one is free to choose 
the rest of this operator to satisfy the constraint © which means that the curve should 

n 

be a closed one. As an explicit example ^7], if one wants to implement a one qubit gate 
u = E U (2) one should take k = 2 and can take n = 3 and X as 



(10) 



where w is a two dimensional vector which should be determined by satisfying the constraint 
(0). In the dimension n is always taken to be A; + 1 and solutions for some one quibt 
gates (e.g. the Hadamard gate) for which k = 2 and two qubit gates (e.g. the CNOT 
gate and the 2 dimensional fourier transform) for which k = 4 are obtained by extensive 
numerical search in the space of two and four dimensional vectors respectively. 
Our aim in this paper is two-fold. First we provide an exact method (free of numerical 
calculations) for obtaining solutions for arbitrary unitary gates acting on an arbitrary 
number of qubits. We obtain in a much simpler method the solutions already obtaind in 
Second we address the problem of scalability of this construction. That is, we modify 




the method of |7[ so that we can holonomically implement any local quantum gate on an 
array of qubits by appropriate choices of curves in the parameter space. This will give us 
full power for doing holonomic computation on a quantum computer in a scalable way. 
The structure of the paper is as follows: In section IH] we present our exact solution for 
the case when n = k + 1. Since this is the case already studied by the authors in 0], who 
have supplied various examples by searching in the space of matrices numerically, we suffice 
to give only proofs for general gates and refrain from giving examples. In section IIIII we 
generalize our method and show how to holonomically implement any local quantum gate 
on an array of qubits by appropriate choices of curves in the parameter space. 
We end the paper by a short discussion. 



II. ANALYTICAL SOLUTION 

Let us take an arbitrary anti-hermitian operator X G u{n). This determines the adjoint 
orbit via equation (jlj). Due to equation (jHJ and the fact that n = k + 1, the general form of 
X is as follows: 



X 



(11) 



where s is a real number and i = Note that the upper k hj k diagonal block of X is 

dictated by the choice of unitary operator U = e""^ that we want to represent holonomically. 
We have choices only for the vector w and the number s. We can use this freedom to solve 
the constraint (jSJ which requires that the loop be closed. 

In the parameter s has been taken equal to zero. However we will see that a better choice 
for this parameter exists which will simplify the final results considerably. Our task is now 
to exponentiate the operator X. 

Exponentiation of a specific anti-hermitian operator is a straightforward (albeit highly te- 
dius) task which can be done by diagonalization of the operator (since we have to solve 
cubic equations with general coefficients). Since we want to do this for a completely general 
operator of the form (|TT|) . for which we have no a priori information about the matrix A and 
the vector w, we make the simplifying restriction that w is proportional to an eigenvector 
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of A, that is 



Aw — iXw, w^w — a 



(12) 



where A is a real parameter. 

With this assumption we can use a recursive method for calculation of e*^ instead of nu- 
merical methods. 

A simple calculation for the first few powers shows that the general form of X" is 



A^ — bnWW^ Cn 



W 



(13) 



where a{n),b{n) and c{n) are coefficients to be determined. These coefficients satisfy the 
following initial conditions: 



60 = 61 = 0, 
Co = ci = 1, 
do — 1 di — is. 



(14) 



To obtain these coefficients we use 



X 



n+l 



W 



^ A^ - bnWW^ CnW ^ 



V 



—W I s 
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from which we obtain the following recursive system of equations 



(15) 



hn+l = i\hn + C* , 
Cn+1 = ^ACn + dn, 
dn+1 isdji Ct Cyi- 



(16) 



To solve this system we define the following generating functions 



n=0 

00 
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n=0 
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(17) 
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Once we find these generating functions we liave tlie exponential as 

' e^^-B{t)ww^ C{t)w 



\^ -C*{t)w^ D{t) 
In terms of generating functions the recursion relations read as follows: 



(18) 



B'{t) = iXB{t) + C*{t), 
C'{t) = iXC{t) + D{t), 

D\t) = isD{t)~a^C{t). (19) 
The second and the third equations can be combined to give 

C"{t) = i{X + s)C'{t) + (As - a^)C(t), (20) 

the solution of which is 

C(t) =7ie''^* + 72e''=*, (21) 

where 

Qi,2 = ^(A + ^ ± V(A-s)2 + 4«2) =: ± u), (22) 



where we have defined u := |y (A — s)^ + Aa^ and 71^2 are constants to be determined from 

initial conditions. 

In view of the initial conditions 

Co = C(0) = and Ci = C"(0) = 1, (23) 

we find the final form of the generating function 

e«i*-e«2t (A+.)t sin i^t 

Cit) = = 2 . (24) 

qi-q2 V 

From this solution we find that a simplifying choice for the parameter s is to take s = —A 
for which we will have: 

m = (25) 



where now 



I, = VA2 + a2 (26) 

. Equation (fT^ will then give 

D{t) = C'{t) - iXC{t) = cos ut - sin ut. (27) 



Finally we solve the first equation of ()19|) by rewriting it as 

-^(e"'^*5(t)) =e-*^*C*(t), (28) 

the solution of which is 

B(t) = f\'^^'-^^C*{T)dT 
Jo 

1 / iXf .A . 

= 1^ 6 — cos z^t — sm z/t 

[i'^ — X^j \ V 

where we have used the initial condition = B{0) = 0. 
What we now require is that the constraint (0) be satisfied. This constraint is satisfied 
when the operator becomes block-diagonal, since in this case [e"^,Po] = 0. Due to 
equation (fTHj) this requires that we set C(l) =0 which in view of demands that we fix 
u as u = Tin or in view of equation ()26|). = (tt-tt)^ — A^, where n is an integer. 
Thus the recepie is as follows: for any given operator U = e~^, we choose an eigenvector v 
of A, we form the matrix X as in |TT| with w = av where a is given as = (nyr)^ — A^. 
Thus for any operator in a A; dimensional space there are a number of k closed curves, each 
corresponding to one of the eigenvectors of A. In view of equation (|24p . for each fixed n 
the Hamiltonian H{t) turns back to its original form Hq exactly n times during the period 
from t = to t = 1. The optimal choice for n is thus n = 1. 




III. GENERALIZATION AND PHYSICAL REALIZATION 

Let us summarize our findings up to now. If we want to implement a unitary operator 
(a quantum gate) on a A; dimensional Hilbert space, we take a system whose Hilbert space 
is A; + 1 dimensional. On this larger Hilbert space, we take a Hamiltonian Hq = Pq, where 
Pq is the projector operator on the smaller k dimensional system and then adiabatically 
transform this initial Hamiltonian by the adjoint action H{t) = e^^ Poe~^^ , where X is a 
vector in the lie algebra G u{k + 1). For any given unitary operator U & U{k), we have a 
suitable vector X, such that the holonomy operator which is adiabatically induced on the 
degenerate k dimensional Hilbert space coincides with that operator, that is e'^c"^ = U, 
where m\A\j{t)) = m\d\j{t)) = {t\X\j), 1 < 2, j < k. 

Actually for any given operator U, we have k generally different solutions X, since we have k 
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different choices for the vector w to insert in , each proportional to a different eigenvector 
of the operator A. This is an important point which we will come to later in this section. 
In fact we will see that it is best to combine the separate solutions into one single solution. 
The reason is that when it comes to physical realization of the technique in section (pi]). 
the problem of dimension becomes an annoying obstacle, since if we want to implement a 
single qubit gate, we require a three dimensional Hilbert space and if we want to implement 
a two-qubit gate, we need a five dimensional Hilbert space and so on. Generally we do not 
have such a control over the dimensions of Hilbert spaces in most experimental scenarios for 
quantum computation. In most of these scenarios we have a number n of two dimensional 
systems (e.g. spins, ions , etc.) which play the role of qubits which when put together 
make a Hilbert space of dimension 2". When one enlarges this system by adding a quibt (an 
ancilla qubit) the dimension of the Hilbert space doubles to 2"+^ instead of getting increased 
by one. 

The most ideal situation is to have a main array of qubits numbered from 1 to n together with 
an ancilla qubit which we number as n + 1. This ancilla qubit is then used for implementing 
holonomically a set of universal gates (one qubit and two qubit gates) on any qubit or pair 
of qubits in the main array of qubits. 

In order to realize this ideal situation we modify our method as follows. We let the dimension 
of the larger Hilbert space to be twice the dimension of the degenerate subspace. For 
simplicity of presentation we first consider the example where we have a single qubit as our 
main system and a single ancilla qubit. Instead of (|TT|l we now take X as follows: 



where wi = aiVi and W2 = OL2V2 are proportional to the two eigenvectors of the operator 
A corresponding to eigenvalues Ai and A2 respectively. In this way we incorporate the two 
separate solutions into one solution for the sake of scalability of holonomic computation. 
We can easily repeat all the calculations in section (^. For example we will have 



A 



{'Wi,W2) 




J 



(30) 



X 



n 




(31) 
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where the functions bk,n, Ck,n and dk,n for = 1, 2 are found to satisfy the folfowing recursion 
relations: 



bk,n+l — i^kbk,n + C^,ni 
dk,n+l = —i^kdk,n — Oilcl ^. 



(32) 



This shows that the recursion relations for various indices k are decoupled and thus the 
solutions for the generating functions found in section (jll} can be directly carried over to 
here. For example we have for each k, that 

sin Ukt 



Ck{t) :-- 



I'k 



(33) 



where Vk = + <^k- This tells us that the matrix will be block diagonal if we choose 
Pk = nkiT or optimally if Vk = tt. 

A more concise and instructive way of doing the above calculation is to rewrite 

{Wi,W2, ■ ■ -Wn) as 

(twi, W2, ■■■Wn) = {aiVi, a2V2, ■ ■ ■ anVn) = ^D, (34) 

where fl is the unitary matrix which diagonalizes A and D = diagonal(o;i, a2, ■ ■ ■ Then 
we can write X in the form 
/ 



X 



V 



A QD 

-Dn^ -iA 




(35) 



In this form the matrix X is easily exponentiated to give 



Q 
/ 



^]t 
/ 



(36) 



where the matrix Xn has been defined as 



Xo-.-- 




(37) 



It is clear that 



Xo = i {a^ (g) A + cTy (g) D) 
Xl = -I®{A^ + D^). 



(38) 



9 



Note that both A and D are diagonal matrices. This then leads to the following result 

Sill i^i 

e^o* = I0cosut + i{a^ A + CTj, ® D){I ® — — ), (39) 



where u := yA^+^D^ is a diagonal matrix which generalizes and replaces the simple factors 
u = vAM-c? of section (jll]). The condition is now satisfied if we choose 



^Xl + al = rifeTT VA;. (40) 



If we take all these factors equal to n which is the optimal choice for each factor, the matrix 
u will equal nl. This brings about further simplifications as we will see. 
From (jnni) and (fH^jl . we find that 



^ cos 7Tt+- sin nt — sin vrt 



V 



■■^^sinTTt cosTTt — i- sinyrt 



(41) 



Now that we have obtained the general form of the operator X, for any unitary operator, 
let us see what is its form for a one qubit gate u := e"' acting on one of the qubits, say qubit 
k, in the main array. Then we will have 

A = I^^-^ ®a® J®"" , n = I^'-' (g)w(g) J^"" , 

A = /^'^-^ ® A ® J®""' , D = I^""-^ ® d ® J®""' , (42) 

where w is the matrix which diagonalizes a, A is the diagonalization of a, and d is the 



diagonal matrix d = diag{ai,a2), where tt = \JX\ + a|. Thus we see that for a local one 
qubit gate uj. = I®^^'^ ® u® J®"~^ which acts trivially on all other qubits except the fc— th 
qubit, our matrix X will also be a local matrix 

Xfc :=/^'=-^(8)z® J®"-^ (43) 

is a matrix acting on the k — th qubit and the ancilla and x is a four by four matrix of the 

a zud 

—dvj^ —iX 

By the same analysis, we find that a local two qubit gate can also be implemented holo- 
nomically, by a local choice of X acting only on the two qubits and the ancilla. The general 
form of the matrix X is always as given in 

For Hq we can assume the following form Hq = |(/ + (Jz) ® I, where the first factor acts 
on the ancilla and the second factor acts on all the other qubits. Thus the Hamiltonian 
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form (jHSjl . that is x 



Hq describes an array of qubits in which only the ancilla qubit is in a magnetic field in 
the z direction. The eigenstates of are |0) (S> \ek),k = 1 to 2" with energy E = 1 
and |1) ® |efc), A; = 1 to 2" with energy E = 0, where \ek) are the basis vectors of the 
main array of qubits. Thus the Hamiltonian is the projection operator on the main array. 
According to ((Tj), one now moves arround the orthonormal frame {|0) |efc)}, by the action 
of the operator e*^ which means that the Hermitian operator iX acts like a Hamiltonian 
in the fixed frame. After one cycle of revolution from t = to t = 1, the desired operator 
corresponding to the choice of X has been implemented holonomically on the main array. 
For a one qubit gate, the operator x can be written in terms of Pauli matrices in the form 
1 i 

X = -(/ + cr^) ® a - -(/ - az) X + (J+ ® zud - a- ® dvj^ . (44) 

Here the first factor of the tensor product acts on the ancilla and the second factor acts on 
a local qubit. As a concrete example if we take u = az, for which a = y(/ — cTz) and zu = I, 
and follow the above construction, after straightforward calculations we find the following 
form of the operator x: 

ITT 

x = — [{(Jz -ay)®I+ {(jy -I)® a,\ . (45) 

IV. DISCUSSION 

We have found exact solutions for the holonomic implementation of local quantum gates 
on an arbitrary array of qubits plus an ancilla. A family of isospectral Hamiltonians all 
with the same spectrum of the projection operator on the main array (without the ancilla) 
is constructed by the choice of a suitable curve in the form of Hit) = e^^ Hqc''^^ , defined 
by an operator in the full Hilbert space of the array and the ancilla. For every single 
qubit gate acting on the A;— th qubit, there is a local operator X^ acting nontrivially only 
on the ancilla and the qubit in question, which implements that gate on the k—th qubit 
holonomically. This is also the case for two qubit gates. 

Thus by imbedding the appropriate operators for a suitable set of one and two qubit gates 
which comprize a universal set of quantum gates, in various places in the tensor product of 
the array, one can enact holonomically any set of gates on any subset of qubits. 
This of course leads us to the actual physical implementation of such a method in a concrete 
physical realization of quantum computers. We leave this problem for future research. 

11 



Acknowledgement We hereby thank A. Sadrolashrafi and A. T. Rezakhani for valuable 
discussions. After submitting this paper we were informed by the authors of [3] that they 
have announced a similar solution (without proof) in their revised paper. 



[1] Zanardi P. and Rasetti M., Phys. Lett. A. 264, 94(1999). 

[2] Ellinas D. and Pachos J., Phys. Rev. A. 64, 022310(2001). 

[3] Pachos J., Zanardi P. Rasetti M., Phys. Rev. A. 61, 010305(R)(1999). 

[4] Pachos J., and Zanardi P., Int. J. Mod. Phys. 15, 1257(2001). 

[5] Pachos J. and Chountasis S., Phys. Rev. A. 62, 052318(2000). 

[6] Fujh K. , Rep. Math. Phys. 48, 75 (2001). 

[7] Tanimura S., Hayashi D., and Nakahara M., Exact solutions of holonomic quantum computa- 
tion, quant-phys/0312079\ v2. Physics Letters A, Vol 325 No 3-4 pp 199-205. 



12 



